Binding and dimerization of PGLa peptides in anionic lipid bilayer studied by replica exchange molecular dynamics

The 21-residue PGLa peptide is well known for antimicrobial activity attributed to its ability to compromize bacterial membranes. Using all-atom explicit solvent replica exchange molecular dynamics with solute tempering, we studied PGLa binding to a model anionic DMPC/DMPG bilayer at the high peptide:lipid ratio that promotes PGLa dimerization (a two peptides per leaflet system). As a reference we used our previous simulations at the low peptide:lipid ratio (a one peptide per leaflet system). We found that the increase in the peptide:lipid ratio suppresses PGLa helical propensity, tilts the bound peptide toward the bilayer hydrophobic core, and forces it deeper into the bilayer. Surprisingly, at the high peptide:lipid ratio PGLa binding induces weaker bilayer thinning, but deeper water permeation. We explain these effects by the cross-correlations between lipid shells surrounding PGLa that leads to a much diminished efflux of DMPC lipids from the peptide proximity at the high peptide:lipid ratio. Consistent with the experimental data the propensity for PGLa dimerization was found to be weak resulting in coexistence of monomers and dimers with distinctive properties. PGLa dimers assemble via apolar criss-cross interface and become partially expelled from the bilayer residing at the bilayer-water boundary. We rationalize their properties by the dimer tendency to preserve favorable electrostatic interactions between lysine and phosphate lipid groups as well as to avoid electrostatic repulsion between lysines in the low dielectric environment of the bilayer core. PGLa homedimer interface is predicted to be distinct from that involved in PGLa-magainin heterodimers.

www.nature.com/scientificreports/peptide:lipid ratio (1:200) PGLa adopts a surface bound S-state with almost parallel orientation of the peptide with respect to the DMPC/DMPG bilayer and upward orientation of lysine side chains 15,16 .However, at higher peptide:lipid ratio (1:50) PGLa shifts to a tilted T-state, in which the C-terminus inserts into the bilayer [15][16][17][18] .In fact, oriented circular dichroism (CD) spectra provided an estimate of the free energy of PGLa homodimers 19 .PGLa is also well known to form heterogeneous aggregates with another AMP magainin resulting in a synergistic increase in antimicrobial action [20][21][22][23] .Although the precise heterodimer structure remains unknown, a combination of experimental techniques suggested that glycine PGLa amino acids Gly7 and Gly11 play a critical role in PGLa-magainin aggregation.It is also of note that PGLa states depend on bilayer curvature 24,25 , thickness of bilayer hydrophobic core 17 , and temperature 26 .Specifically, an increase in temperature induces the PGLa transition from T to S states in the liquid crystalline phase.Binding of PGLa to the DMPC/DMPG bilayer has been shown to lead to significant disordering in lipid structure 27 .
In parallel to experimental investigations, molecular dynamics (MD) simulations have been used to probe PGLa interactions with lipid bilayers.Recent 1 μs long MD simulations based on AMBER force field have revealed the helical secondary structure and orientation, including tilt and side chain rotation angles, of PGLa monomers in the DMPC bilayer to be in good agreement with the NMR analysis 28 .Separate AMBER MD simulations showed that the orientation of L18W-PGLa peptide in POPE/POPG bilayer is consistent with X-ray and neutron scattering data 21 .It followed from the MD simulations utilizing CHARMM force field that PGLa monomers bound to the DMPC membrane maintain helical structure and reside just below the phosphorus lipid atoms 29,30 .The analysis of preconstructed antiparallel PGLa dimers has indicated that such assemblies have higher stability than parallel dimers 29 .In order to probe the conformational ensemble of PGLa peptides, our group has performed replica exchange with solute tempering (REST) simulations of PGLa binding to anionic DMPC/DMPG bilayers 31,32 at the low peptide:lipid (P:L) ratio.By analyzing the binding free energy landscape we identified two major bound states, a metastable surface bound state and a dominant inserted state.In both states positively charged PGLa amino acids maintain electrostatic interactions with anionic phosphate groups by rotating the PGLa helix around its axis.PGLa binding causes an influx of anionic DMPG and efflux of zwitterionic DMPC lipids from the peptide proximity.The analysis of binding free energy suggested that PGLa binding to the DMPC/DMPG bilayer is governed by the balance between desolvation of PGLa positive charges and formation of electrostatic PGLa-lipid interactions.
In this paper, we extended our REST investigation to the PGLa peptides binding to the DMPC/DMPG bilayer at the high P:L ratio.Computation of free energy landscapes suggests that PGLa homodimers are diverse and metastable existing in dynamic equilibrium with monomeric species.Importantly, the PGLa monomers occurring at the high P:L ratio are distinct from those observed at the low ratio.Unexpectedly, PGLa interpeptide interactions reduces the peptide helical propensity and forces the dimerized peptides to position themselves closer to the bilayer surface adopting a criss-cross dimer interface.We also found that the high P:L ratio results in profoundly weaker efflux of zwitterionic lipids from the PGLa binding footprint.The physicochemical mechanism, which changes the properties of PGLa at the high P:L ratio, is proposed.

PGLa properties and bilayer structure at the high peptide:lipid ratio
PGLa properties within the DMPG/DMPC bilayer: Our REST simulations examined PGLa peptides binding to the anionic DMPC/DMPG (3:2) bilayer at the peptide:lipid (P:L) ratio of ∼ 1:25 (Fig. 1).At 330 K PGLa binds to the bilayer and dimerizes with the probabilities P b ∼ 1.0 and P d = 0.29 ± 0.05 , respectively.Therefore, the simulations produced the conformational ensemble composed of a mixture of bound PGLa monomers and dimers.Because multiple PGLa species, monomers or dimers, coexist, we refer to them collectively as those sampled at the high P:L ratio.As a reference we use the PGLa binding to the same bilayer at 310 K and lower P:L ratio 1:50 reported previously 32 .The use of "high" and "low" P:L ratio terms for the simulations requires caution.Due to simulation design at the low P:L ratio PGLa cannot aggregate, whereas at the high P:L ratio the peptide aggregation is limited to dimers.Therefore, an alternative reference to the high and low P:L ratio simulations could be 2-or 1-peptide systems, respectively.We discuss the correspondence between the experimental and computational P:L ratios later in this section.Also, in Supplementary Information (SI) we show that the temperature dependence of PGLa properties in CHARMM force field is weak.This circumstance allows us to compare the PGLa properties observed at the high P:L ratio and 330K with those at the low P:L ratio and 310K.
Figure 2 shows that at the high P:L ratio the helical propensity H(i) is reduced for all PGLa amino acids.Indeed, the overall helical fraction H is 0.40 ± 0.06 at the high P:L ratio compared to 0.63 ± 0.05 at the low, i.e., it is reduced by more than a third.Furthermore, at the low P:L ratio 17 PGLa amino acids are classified as helical for their H(i) > 0.5 , whereas at the high ratio their number decreases three-fold to 6.The helical pro- pensity is particularly suppressed in the C-terminus, where H(Ct) reaches 0.81 ± 0.05 at the low P:L ratio, but drops to 0.51 ± 0.07 at the high ratio.Besides helix at the high P:L ratio, PGLa forms turns and random coil with their fractions being 0.49 ± 0.06 and 0.10 ± 0.02 .Multiple experiments have reported the formation of helix in PGLa.According to NMR data probing PGLa binding to 3:1 POPC/POPG micelles at the low 1:100 P:L ratio, the peptide is mostly helical from the position i = 6 to 21 12 .Similarly, CD studies showed that upon binding to DMPC/DMPG bilayer the PGLa helical fraction H varies from 0.67 to 0.76 in the temperature range from 303 to 318K at the P:L ratio of 1:200 10 .A more recent CD study of a similar system at the P:L ratios of 1:200 and 1:50 reported H to be about 0.6 14 .These findings are in good agreement with our simulations of PGLa monomer 32 .Raman spectroscopy experiments reported PGLa secondary structure upon binding to DMPC liposomes at much higher P:L ratio of about 1:3 and ambient temperature 33 .In contrast to the data at lower P:L ratios, the PGLa helical structure propensity was strikigly lower reduced to 0.42 ± 0.05 .In qualitative agreement, our simulations report that at the elevated P:L ratio < H > dropped to 0.40 ± 0.06 .There are no experimental data on the helical fraction at the P:L ratio exceeding 1:50.Thus, it appears that the PGLa helix structure is sensitive to P:L ratio and partially unravels with its increase.We have developed a rationale for decreasing helical fraction with the increase in P:L ratio, which is proposed after the analysis of lipid distributions.Because it is challenging to unambiguously match the experimental and simulation P:L ratios (see below), our prediction that < H > decreases with the P:L ratio should be viewed as preliminary, and our simulations may still underestimate < H >.
To map the PGLa location within the DMPC/DMPG bilayer, we first computed the probability distribution P(z CM ) of the PGLa peptide center of mass z CM along the bilayer normal.Figure 3 reveals that with the P:L ratio increase the peptide shifts deeper into the bilayer.In fact, the average position of the PGLa center of mass z CM changes from 13.1 ± 0.9 at the low to 11.6 ± 0.5 Å at the high P:L ratio.To get further insight, we used the probabilities P(z; i) for amino acid i to occur at the distance z from the bilayer midplane to compute the average positions of amino acids z(i) (see Methods).We found that with the increase in P:L ratio there is a discernible deeper insertion of PGLa C-terminal amino acids, and at both P:L ratios all amino acids, on an average, reside below the position of the center of mass of phosphorus atoms in a leaflet z P .Thus, at both P:L ratios all PGLa amino acids are classified as inserted (see Methods).Previous microsecond-long simulations, which did not utilize enhanced sampling, also found that PGLa monomers are deeply inserted in the DMPC bilayer residing well below phosphate groups with their center of mass z CM = 9.4 Å 29 .
The orientation of the PGLa peptide in the bilayer has been studied experimentally using the tilt angle γ (see Methods).Because the N-or C-termini exhibit a drop in the helical propensity (Fig. 2), we computed γ defined for the peptide region 6-14.(The results do not appreciably change if γ is computed for the region 15-20).Fig- ure 4a compares the probability distributions P(γ ) computed at the high and low P:L ratios.Both distributions are approximately unimodal with the average �γ � = 107 ± 5 • and 89 ± 2 • , respectively.Next, we considered the rotation of the aliphatic face of the PGLa helix.To this end, we used the rotation angle β (see Methods) defined for Lys19 and computed the probability distribution P(β) . (Qualitatively similar results have been collected if β is defined for Lys15 or Lys12).Figure 4b shows that at the high P:L ratio Lys19 side chain predominantly points up with the average angle �β� = 105 ± 5 • , and since the peptides are mostly inserted, it is directed up toward the anionic phosphate groups.At the low P:L ratio β is 74 ± 12 • .Thus, changing P:L ratio increases the peptide tilt and maintains Lys19 upward direction.The experimental values for the tilt and Lys12 rotation angles have been  reported for the mixed DMPC/DMPG bilayer at 308 K, full hydration, and different P:L ratios 15,16 .At the 1:200 ratio the experimental γ = 98 • and β = 115 • .However, when the P:L ratio increases to 1:20 γ and β become 121 • and 113 • .
Correspondence between simulation and experimental conditions: To properly compare the simulations and experiments, we need to find the correspondence between the experimental and simulation P:L ratios.Our low P:L ratio simulations, which keep PGLa exclusively monomeric, may roughly correspond to the experimental P:L ratio of 1:200, because our helical fraction < H >= 0.63 and the PGLa tilt angle < γ >= 89 • agree approxi- mately with the experimental H = 0.60 and γ = 98 •14,15 .It is more challenging to match our high P:L ratio to the experiment, because by design our simulations do not allow aggregation beyond a dimer.To this end, we computed the quadrupolar splittings �ν q (see SI). Based on the direct comparison of in silico and experimental �ν q it appears that our high P:L ratio represents the experimental 1:20 P:L ratio (see Table S1) 15 .Incidently, this outcome matches the actual P:L ratio (1:25) in our simulation system.However, some difference in the tilt γ and rotation β angles between our simulations and experiments must still be noted (see above).These discrepan- cies can be explained by the differences in angle definitions, difficulties in precise matching experimental and computational conditions, and differences in bilayer composition.Nevertheless, experiments and simulations are in qualitative agreement with respect to data trends, namely, that with the increase in P:L ratio the PGLa helical propensity decreases, while the overall PGLa tilt increases.Simulteneously, the upward direction of Lys side chains in the PGLa helix is preserved.
PGLa impact on the DMPC/DMPG bilayer: The changes in the DMPC/DMPG bilayer structure caused by PGLa binding were evaluated using the volume number density n l (r, z) of lipid heavy atoms computed as a func- tion of the distance r to the peptide center of mass and the distance z to the bilayer midplane (see Methods).The probability distributions P(β) describe the orientation of PGLa helix aliphatic face.Data in black and grey in both panels correspond to the high and low 32 P:L ratios.This figure shows that PGLa tilt increases at the high P:L ratio.
Figure 5a compares n l (r, z) for the high and low P:L ratios.In both cases, PGLa binding creates a lipid density void around the peptide, but its extent depends on the P:L ratio.If the bilayer thinning D is defined as the difference in the average bilayer thickness D in the near and distant regions, then it is 6.3 ± 0.9 and 17.3 ± 0.3 Å at the high and low P:L ratios, respectively.The surprising result is that PGLa binding at the high P:L ratio causes much weaker bilayer thinning than at the low.A possible rationale is proposed below after the analysis of surface lipid density.
Our previous study at the low P:L ratio has demonstrated that PGLa induces a water invasion into the DMPC/ DMPG bilayer 32 .To check this effect at the high P:L ratio, we computed the probability P ww (z) for a water wire to penetrate the bilayer to the depth z or below (see Methods).Figure 5b compares P ww (z) obtained at the high and low P:L ratios.The figure unambiguously show that at the high P:L ratio water invades the bilayer deeper.
In fact, at the low P:L ratio half of water wires extends at least to the depth of z = 11.5 Å compared to z = 8.5 Å at the high P:L ratio.Thus, the increase in the P:L ratio facilitates water penetration in the bilayer by additional 3 Å.This outcome can be readily understood if we recall that the average position of the PGLa peptide center of mass z CM at the high P:L ratio is 2.5 Å deeper than at the low (Fig. 3a).

(a) (b)
High P:L ratio Low P:L ratio Continuous black lines mark the bilayer boundaries z b (r) .The boundaries of near, proximal, and distant regions are given by dashed lines.The density distribution n l (r, z) is computed by permuting the upper and lower leaflets and peptides placing the reference peptide center of mass at r = 0 .This figure implicates surprisingly weak disruption of the bilayer structure at the high P:L ratio.(b) The probabilities P ww (z) map the insertion of the water wires to the depth z or below.Data for the high and low P:L ratios are in black and grey.
In our previous studies at the low P:L ratio, we observed strong interactions between cationic PGLa and anionic DMPG lipids resulting in an influx of these lipids into the peptide proximity with simultaneous efflux of zwitterionic DMPC 31,32 .To examine the same effect at the high P:L ratio, we considered the lipid surface number densities n s,x presented in Table 1.At the high P:L ratio PGLa decreases n s,l in the near region compared to the distant one by merely 16% compared to the two-fold drop at the low ratio.Decomposing densities in Table 1 reveals a negligible increase in the near DMPG surface number density n s,DMPG and a 37% drop in the DMPC density n s,DMPC .At the low P:L ratio there is also a minor variation in n s,DMPG , but, in contrast, a very sharp, three-fold drop in n s,DMPC .Furthermore, at the high P:L ratio the DMPC:DMPG ratio in the near region is 1.2, whereas it is 0.5 at the low P:L.Hence, at the high P:L ratio the efflux of DMPC lipids from the near region, although present, is nevertheless much muted than at the low ratio.In fact, at the high P:L ratio the concentration adjusted DMPG density is still about one-third less than n s,DMPC , whereas it surpasses n s,DMPC at the low ratio (Table 1).The plausible reason for this effect is as follows.Table 1 shows that at the low P:L ratio the proximal region exhibits a minor drop in DMPC density compared to the distant region.Then, if at the high ratio the average interpeptide distance is 14.8 ± 0.5 Å, a DMPC lipid occurring in the proximal region of one peptide becomes a "near" lipid for the other.The outcome is a cross-correlation between peptide lipid shells that partially levels off the densities in the near and proximal regions.
The analysis of lipid densities provides an insight into the bilayer thinning.Because at the high P:L ratio the correlations between PGLa lipid shells largely sustain the lipid density near and away from the bound peptides, the bilayer thinning is muted compared to the low P:L ratio.In fact, according to Table 1 the near density n s,l at the high P:L ratio is increased by about 50% compared to that at the low ratio.Furthermore, a weak efflux of DMPC lipids offers a rationale for the reduced helical fraction in PGLa peptides at the high P:L ratio.Because the near bilayer region becomes less anionic, it induces weaker helical propensity in the PGLa peptide due to smaller polarity gradient between the bilayer surface and interior.Partial unraveling of the helix allows lipids to enter the PGLa volume reducing the overall drop in the near lipid density and further alleviating the bilayer thinning.We are not aware of experimental studies probing these effects, apart from the report that PGLa-induced clustering of anionic lipids has been observed in differential scanning calorimetry experiments 34 .Notably, no clustering of DMPG lipids has been detected around another AMP, indolicidin 35 .Although its smaller net positive charge compared to that of PGLa may offer a straightforward reason, another, more subtle contribution cannot be ruled out.Indolicidin experiments were performed at the P:L ratio of up to 1:10, whereas PGLa experiments have used 1:20.Our data argue that high P:L ratio may impede DMPG clustering.

Free energy landscape of PGLa binding and dimer formation
According to our previous studies at the low P:L ratio PGLa monomers bind to the DMPC/DMPG bilayer via a two-state mechanism.Specifically, the PGLa peptide adopts a dominant inserted state and a metastable surface bound state 32 .Then, to investigate the free energy of PGLa binding and dimer formation, we selected two reaction coordinates -the position of the center of mass of two PGLa peptides in a leaflet z DCM and the number of inter- peptide contacts C d (see Methods).Figure 6a presents the free energy of PGLa peptides in a leaflet G(z DCM , C d ) computed as a function of these two variables.The free energy landscape reveals a thermodynamically stable state, inserted monomers (IM), augmented by a manifold of metastable dimeric states.In Fig. 6a one can distinguish the nascent dimers, which are not separated by any free energy barrier from IM and represent transient assemblies readily falling apart and converting back to IM.To partition nascent from "mature" dimers D, we identified the set of saddle points or transition states on the free energy surface G(z DCM , C d ) in Fig. 6a.Consequently, we focus our analysis on IM and D, whose characteristics are listed in Table 2.
According to Table 2 the most thermodynamically stable state, the monomeric IM, occurs with the probability 0.71 ± 0.05 and is displayed in Fig. 6b.It follows from Fig. 7a that IM peptides are deeply inserted into the bilayer hydrophobic core positioned, on average, at z CM = 10.6 Å from the midplane.IM helical fraction is 0.42 and according to Fig. 7b stable helix appears at most positions i ≥ 6 .According to Table 2 the peptide is weakly tilted toward the bilayer core, although Fig. 7c suggests a broad distribution of tilt angles with the standard deviation δγ = 19 • .Since IM rotation angle β in Fig. 7d exhibits a unimodal distribution around the average of 103 • , the helix approximately positions its cationic Lys19 side chain upward in the direction of anionic phosphate groups.It is important to note that IM differ from the monomeric states observed at the low P:L ratio.www.nature.com/scientificreports/originates from the cross-correlation between lipid proximal and near regions surrounding the peptides at the high P:L ratio as described above.Since IM species are dominant, the increase in PGLa tilt at the high P:L ratio is primarely due to IM.Similarly, the weak bilayer thinning at the high P:L ratio is caused by IM species and better mixing of lipids and amino acids at the high P:L ratio.
The dimers D occur with the probability 0.15 ± 0.04 and are separated from IM by the free energy gap of 1.0 kcal/mol.Importantly, a high free energy barrier of almost 4 kcal/mol occur on the IM→D path, whereas a much lower barrier exists along D→IM (Table 2).Thus, PGLa D dimers once formed can be dissociated into IM monomers.The free energy G(z DCM , C d ) in Fig. 6a indicates that D are diverse species with broad distribu- tion of insertion depths z DCM and the number of interpeptide contacts C d .The representative D dimer structure is shown in Fig. 6c.The average position of PGLa center of mass in D z CM is 16.0Å, which is more than 5Å above the average location of IM peptides (Table 2).Indeed, compared to IM the dimer probability distribution P(z CM ) in Fig. 7a exhibits a striking shift to the bilayer surface.Thus, in contrast to the monomeric species the D dimers are partly expelled from the bilayer.It follows from Fig. 7b and Table 2 that apart from a moderate destabilization of helix at the positions 6 through 11 the helical propensity in D is generally similar to that of Table 2. PGLa peptide states.a Fraction of PGLa peptides in a state k.To find P(k), we included all structures associated with k independent of their The transition k→l crosses the free energy barrier IM.According to Fig. 7c the probability distribution P(γ ) for D is broad, but it is still apparent that the D pep- tides have smaller tilt than IM.Indeed, according to Table 2 the average D tilt γ is 93 • whereas it is 110 • for IM.Furthermore, the probability distributions of rotation angles P(β) in Fig. 7d and the average β for D and IM are also distinct.Although D P(β) is broad, it is extended toward larger values compared to IM.As a result Lys19 side chain direction deviates, on an average, by almost 40 • (50 • for Lys12) from the bilayer normal, whereas in IM the respective deviation is merely 13 • .In summary, there are three differences between IM and D. First, D dimers are partially expelled from the DMPC/DMPG bilayer residing near the lipid phosphate groups.Second, the D peptides tend to align parallel to the bilayer surface.Third, D Lys side chains are slanted away from the bilayer normal.It is of note that because of low stability of D the simulation quadrupolar splittings represent the mixture of IM and D states.Next, we examine the dimer interface.The average number of interpeptide contacts in D is 12.8, the most probable of which are listed in Table S2.It is noteworthy that out of 18 contacts occurring with the probability P c > 0.1 , only one (and the least stable) involve Lys amino acid, whereas all others are established between hydrophobic amino acids.In fact, in descending order PGLa dimer formation most frequently involves Ala17, Ala20, Ile13, and Leu21.Furthermore, out of 18 contacts 13 or 72% occur between C-terminal amino acids, but most contacts are off-registry.To determine the D dimer interface geometry, we computed the angle between the peptides α (see Methods).The probability distribution P(α) in Fig. S5 shows a sharp peak at α ≃ 100 • resulting in the average value of 93 • .Thus, although P(α) in Fig. S5 implicates a multitude of peptide alignments in D, overall the dimers prefer an approximately criss-crossed off-registry interface.
The analysis of dimer interface shed light on the differences between IM and D. Because PGLa dimer formation almost exclusively involves apolar amino acids positioned on the opposite face of PGLa helix from the Lys www.nature.com/scientificreports/side chains, the latter ones tend to point outward in the dimer.Consequently, they could not be directed toward the phosphate groups in both peptides if inserted deeply into the bilayer.Our previous studies have implicated the electrostatic interactions between Lys and lipid phosphate groups as a major factor in PGLa binding 31,32 .Then, to maximize the Lys interactions with lipid phosphate groups while preserving the dimer interface, D must position itself close to the bilayer surface.Another reason for the dimer expulsion is that dimer assembly in the bilayer hydrophobic core brings Lys close in a low dielectric environment.Therefore, to minimize electrostatic repulsion while preserving the dimer, it should be repositioned closer to water, which provides better electrostatic screening.The PGLa peptides in D state comply with these scenarios by shifting closer to the bilayer surface and slanting Lys side chains away from the bilayer normal.
Previous simulations have investigated PGLa dimers.In particular, Ulmschneider et al. have performed constant temperature simulations probing the PGLa monomers and dimers 29 .The study reported an increase in the PGLa dimer tilt angle γ to 104-107 • , which is almost identical to our γ .However, upon closer examination this agreement is entirely attributed to the dominant monomeric state IM in Table 2, whereas our PGLa dimers D exhibit, on an average, no tilt.It is of note that our simulations study de novo PGLa dimer assembly using enhanced sampling, whereas previous studies examined the stability of preconstructed dimers.Furthermore, previous simulations suggested that antiparallel dimer interface is the most stable supported by hydrophobic contacts between Leu, Ala, Gly side chains.In contrast, we observed criss-crossed dimer interface geometry and saw no significant contribution from Gly to the dimer interface.This outcome is interesting, because multiple experimental methods pointed to the critical contribution of Gly7, Lys15, and Lys19 to PGLa-magainin heterodimer assembly 20 .A similar conclusion on the importance of salt bridges has been made in the small angle X-ray and neutron scattering experiments 21 .Therefore, our study, which did not find these interactions to contribute to dimer formation, implicates very different interfaces in PGLa homodimer and PGLa-magainin heterodimer.Finally, we can compare the free energies of dimer formation extracted from our simulations and oriented circular dichroism experiments 19 .The probability of PGLa monomers P m in our simulations is 0.71 ± 0.05 , whereas the probability of forming dimers, including mature D and nascent in Fig. 6a, is P d = 1 − P m = 0.29 ± 0.05 .Then, assuming a simple monomer-dimer equilibrium, one can compute the free energy of dimer formation as where V = 59 Åx 59 Åx 37 Å is the simulation volume available for the pair of peptides and c 0 = 1 M is the standard concentration.This estimate is consistent with the experimental value of −2.4 kcal/mol reported for the DMPC bilayer 19 .The agreement should be considered approximate, because our study used anionic DMPC/ DMPG bilayer, and Eq.(1) strictly holds only for large simulation systems 36 .Nevertheless, both simulations and experiments agree that the dimerization propensity of PGLa peptides is weak, and association of PGLa with magainin is required for synergistic antimicrobial activity.Poor optimization of PGLa for forming homodimers may indicate that this peptide perform additional biological functions beyond that of AMP.

Simulation system
The simulation system consisted of the DMPC/DMPG (3:2) lipid bilayer formed by 38 anionic DMPG and 60 zwitterionic DMPC lipids, four PGLa peptides, two of which were placed on either bilayer side, and 4410 water molecules (Fig. 1).The system net charge was set to zero by adding 38 sodium ions and 20 chlorine ions.The dimensions of the system were approximately 59 Å x 59 Å x 74 Å.Thus, the P:L ratio is about 1:25.The DMPC/ DMPG bilayer was selected because it is well studied experimentally 37 and frequently used in the PGLa experiments.The simulation system is symmetric, in which two PGLa monomers bind and dimerize in either bilayer leaflet effectively doubling the conformational sampling.Such design also reduces the pressure disparity in the leaflets 38 .PGLa peptides were modeled using the all-atom CHARMM22 force field with CMAP corrections 39 , whereas the all-atom CHARMM36 force field was used for lipids 40 (see SI).A modified TIP3P model was utilized for water 41,42 .To mimic experiments, the PGLa C-terminus was amidated, while the N-terminus was protonated and positively charged resulting in the PGLa net charge of +5 7 .

Replica exchange simulations
Following our previous studies we used isobaric-isothermal replica exchange molecular dynamics with solute tempering (REST) to probe binding of PGLa to the lipid bilayer 43 .All details of REST implementation can be found in our previous publication 32 .In brief, we considered R = 16 system replicas, which were simulated at the temperatures distributed exponentially from T 0 = 330 to T R−1 = 450 K. REST interaction scaling effectively divides the system into tempered ("hot") and untempered ("cold") partitions.The PGLa peptides along with chloride counterions constituted the "hot" solute, whereas the DMPC/DMPG bilayer, water, and remaining ions were "cold".By restricting tempering to the solute, REST considerably reduces the computational load compared to traditional replica exchange.However, "cold" solute may slow down system equilibration requiring to extend REST simulations or elevate sampling temperature 32 .A possible solution suggested recently by us is to use replica exchange with hybrid tempering (REHT) 32 .However, because REHT requires an increase in R, we opted to use REST, but elevated sampling temperature to compensate for REST deficiencies.Replica exchanges were attempted every 2 ps.
NAMD with REST implementation 44 along with in-house scripts were used to manage the REST simulations.We used periodic boundaries and 1 fs integration timestep.Hydrogen-associated covalent bonds were constrained by SHAKE.To compute electrostatic interactions we applied Ewald summation, whereas force switching was utilized to turn off van der Waals interactions in the interval from 8 to 12 Å.Temperature was controlled by underdamped Langevin dynamics with a damping coefficient γ = 5 ps −1 .To maintain pressure of

Figure 1 .
Figure 1.(a) PGLa peptide sequence.N-and C-terminal regions are shown in pale green and orange.Positively charged amino acids are in blue.(b) Structures of DMPC or DMPG lipids, which are composed of a choline or glycerol group (L1), phosphate group (L2) with phosphorus P atom, glycero backbone (L3), and two fatty acid tails (L4 and L5).L1-L3 represent polar headgroups, whereas L4 and L5 constitute the hydrophobic core.(c) A snapshot of REST simulation system.A pair of PGLa peptides bind to each bilayer leaflet.The bilayer is composed of DMPC/DMPG lipids in pale blue and red, respectively.Peptides are colored following (a).Water is shown as dots, sodium and chlorine ions are purple and yellow spheres, Lys amino acids are in dark blue, and phosphorus P atoms are in tan.PGLa dimer is formed in the upper leaflet.

Figure 2 .
Figure 2.The helical H(i) fractions are computed for PGLa amino acids i. Data in black and grey represent the high and low 32 P:L ratios.Sampling errors are given by vertical bars.The amino acids from Nt and Ct regions are colored following Fig.1a.The figure shows that the high P:L ratio suppresses helical structure.

Figure 3 .
Figure 3.The probability distributions P(z CM ) report the location of the PGLa center of mass z CM along the bilayer normal.Data in black and grey correspond to the high and low 32 P:L ratios.The figure reveals deeper insertion of PGLa peptides at the high P:L ratio.

Figure 4 .
Figure 4. (a) The probability distributions P(γ ) probe the PGLa peptide tilt within the lipid bilayer. (b)The probability distributions P(β) describe the orientation of PGLa helix aliphatic face.Data in black and grey in both panels correspond to the high and low 32 P:L ratios.This figure shows that PGLa tilt increases at the high P:L ratio.

Figure 5 .
Figure 5. (a)The volume number density of lipid heavy atoms n l (r, z) as a function of the distance r to the PGLa center of mass and the distance z to the bilayer midplane maps DMPC/DMPG bilayer structure around bound PGLa peptide.Left and right panels compare the density profiles at the high and low 32 P:L ratios.Continuous black lines mark the bilayer boundaries z b (r) .The boundaries of near, proximal, and distant regions are given by dashed lines.The density distribution n l (r, z) is computed by permuting the upper and lower leaflets and peptides placing the reference peptide center of mass at r = 0 .This figure implicates surprisingly weak disruption of the bilayer structure at the high P:L ratio.(b) The probabilities P ww (z) map the insertion of the water wires to the depth z or below.Data for the high and low P:L ratios are in black and grey.

Figure 6 .
Figure 6.(a) The free energy of PGLa peptides in a leaflet G(z DCM , C d ) computed as a function of the position of their center of mass z DCM along the bilayer normal and the number of interpeptide contacts C d .The contour lines are incremented by 0.5 kcal/mol.The locations of monomeric IM and dimeric D species are shown.The dashed white line approximately partitions nascent dimers and "mature" dimers D. (b) Representative monomer structures are taken from the state IM.(c) Representative dimer structure is taken from the D free energy minimum at z DCM ∼ 17 Å and C d ∼ 16 .The D dimer is positioned at the bilayer-water interface and exhibits a criss-cross interpeptide interface.

SBFigure 7 .
Figure 7.Comparison of PGLa monomers and dimers at the high P:L ratio: (a) The probability distributions P(z CM ) of PGLa peptide center of mass position z CM computed for the monomers IM and dimers D. (b) The helical propensities H(i) for amino acids i are computed for PGLa monomers IM and dimers D. The amino acids from Nt and Ct regions are colored following Fig. 1a.(c) The probability distributions P(γ ) probe the PGLa peptide tilt within the lipid bilayer.(d) The probability distributions P(β) describe the orientation of PGLa helix aliphatic face.In all panels data in blue and red represent monomers IM and dimers D.

Table 1 .
The reason Lipid densities in bilayer regions.